Raman spectra of bilayer graphene to probe the electrostatic environment 
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The Raman shift, broadening, and relative Raman intensities of bilayer graphene are computed as 
functions of the electron concentration. We include dynamic effects for the phonon frequencies and 
we consider the gap induced in the band structure of bilayer graphene by an external electric field. 
We show that from the analysis of the Raman spectra of gated bilayer graphene it is possible to 
quantitatively identify the amount of charges coming from the atmosphere and from the substrate. 
These findings suggest that Raman spectroscopy of bilayer graphene can be used to characterize the 
electrostatic environment of few-layers graphene. 

PACS numbers: 78.20.Bh,63.20.kd,78.30.Na,63.22.Np,81.05.Uw 



I. INTRODUCTION 

Graphcnc-bascd systems have recently attracted much 
attention from both the experimental and the theoreti- 
cal point of view. Graphene is in fact characterized by a 
high carrier-mobility^ii^ which make those systems very 
exciting in view of future applications in nanoelectron- 
ics. In standard experimental setups, few-layers graphene 
(FLG) interacts with the environment through doping 
charges coming from the top (nt op ) and from the bot- 
tom (riboi) of the system. These charges determine the 
external electric field [E = (n top — rtf, ot )|e|/(2eo)] and the 
total electron concentration (n = ntop +^6ot)- E and n 
can thus be varied independently changing the charges 
from the two sides. n top and i%bot can be intention- 
ally induced by applying gate voltage differences between 
the system and the substrate, or via atoms/molecules 
deposition on top of the system. On the other hand, 
important unintentional doping charges are typically 
present. For instance, in FLG obtained by exfoliation 
on Si0 2 ii2i3ASigiL»i'»i"' 1 "i"43 1 i4 1 i5 1 i6 1 i7,ia 1 i9 or ep itaxi- 

ally grown on SiC^ a charge transfer occurs between 
the substrate and the system, giving rise to a finite ribot- 
In analogy, an additional n top can be accidentally induced 
by the uncontrolled adsorption of molecules from the at- 
mosphere. 

Among FLG, the bilayer graphene is particularly in- 
teresting because it becomes a tunable band-gap semi- 
conductor under the application of an electric field per- 
pendicular to the system.— ^ i 20 ' 21 i 22 ' 23 i 24 In this context, 
the determination of the electric field is crucial in order 
to control the band-gap. Moreover, in graphene charge 
impurities originating from the top or from the substrate 
are the main source of scattering which reduces conduc- 
tion performances^ Therefore, the determination of the 
charge transfer from the substrate or from the atmo- 
sphere is highly desirable for possible applications. Al- 
though experimental techniques allow to estimate the to- 
tal electron concentration on the system, an accurate de- 
termination of the respective values of top and bottom 
charges has never been achieved so far, and can be partic- 
ularly challenging, for instance, when doping is intention- 



ally induced by deposition of molecules or polymeric elec- 
trolyte. In this work we use a tight binding (TB) model 
fitted on ab initio calculations to compute the Raman 
shift, broadening, and relative Raman intensity of the G 
modes in bilayer graphene, as a function of the electron 
concentration, for different values of top charges. In par- 
ticular, the screening properties of the system in presence 
of an external electric field are described using ab initio 
density functional theory calculations (DFT), including 
the GW correction, while a TB model is used to repro- 
duce the DFT calculated, GW corrected, band structure 
in the full Brillouin zone. We show that from the mea- 
sured Raman spectra of bilayer graphene it is possible 
to determine the external charge distribution and thus 
the external electric field and the actual carriers concen- 
tration. This result is especially relevant since it shows 
that Raman spectroscopy, which is already widely used 
to investigate graphene-based systems, can be used to 
characterize the electrostatic environment of the sample. 
Moreover, since the charges coming from the atmosphere 
and the substrate are not expected to depend on the num- 
ber of layers, Raman measurements on bilayer graphene 
can also be used to determine the origin and the amount 
of the unintentional doping of monolayer and few-layers 
graphene in the same environment. 

The Raman spectra of monolayer graphene are char- 
acterized by a doubly-degenerate G peak (Eig mode) 
at around 1580 cm _1 i2i2iMiffi This mode shows a 
strong clcctron-phonon coupling, which induces a phonon 
renormalization when n is varied. Therefore, Raman can 
be used to measure the total electron concentration. In 
bilayer graphene, in the absence of an external electric 
field the G peak splits, as in graphite, in a doubly- 
degenerate Raman active mode (-E^g) and a doubly- 
degenerate, Raman inactive, infra-red active mode (E u ). 
The E2 g mode is characterized by a symmetric in-phase 
displacement of the atoms in the two layers (Fig[T]-a), 
whereas E u is characterized by an antisymmetric out-of- 
phase displacement of those atoms (Figfljb). Most Ra- 
man measurements on bilayer graphene show a single G 
peak whose position is used, as in monolayer graphene, 
to measure the electron concentration on the system i 1:L i 14 
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FIG. 1: (Color online) Schematic representation of the Ra- 
man active mode Ei g (a) and the Raman inactive mode E u 
(b) in bilayer graphene. A.B are the inequivalent carbon 
atoms in the top layer, A' and B' are the inequivalent carbon 
atoms in the bottom layer. In a Bernal stacking configuration, 
A and A' are the superposed atoms. 



Interestingly, the splitting of the G peak has been ob- 
served in the Raman spectra of bilayer graphene^ and it 
has been recently attributed to symmetry breaking due to 
the application of an external electric fields Moreover, 
other experimental works recently reported on the infra- 
red spectra of gated bilayer graphene i 5 ' 16 ' 17 ' 18 ' 19 These 
findings suggest that a deep understanding of the behav- 
ior of a splitted G mode in gated bilayer graphene could 
lead to quantify not only the total electron concentration 
n but also the separate values of nt op and n\, t- 



II. COMPUTATIONAL DETAILS 

In order to calculate the electron-phonon coupling 
(EPC) component of the phonon frequencies, broaden- 
ings, and relative Raman intensities as functions of n 
and ntop, we consider the V phonon self-energy^^ pro- 
jected onto the subspacc of the two E 2g and the two E u 
modes: 
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where the sum is on the electron wave vector k and 
the electronic 7r bands ^,^=1,4 are the phonon in- 
dexes, Nk is the number of k vectors, and /ki is the 
occupation of the electron state |ki) with energy e ki . 
D£ = (ki\AH^\kj) is the EPC and AH^ is the Hamilto- 
nian derivative with respect to the atomic displacement 
corresponding to the a phonon. n equals 0.009 cV and 
M is the atomic mass. 

From the phonon self-energy one can calculate only 
frequency variations due to changes in the electron con- 
centration and in the band structure. In order to obtain 



the absolute frequencies we use the following 4x4 matrix: 
n M „(ratop,re) = (u) + Aw(n) - n + i-y^ 1 ) <W 

), (2) 

where wo=1581.5 cm' 1 is the experimentally mea- 
sured frequency of the Raman active G peak in bilayer 
graphene, in absence of doping and electric fieldJ^ 7 an 
is the contribution to the broadening of the G modes 
in graphene and graphite from the anharmonic phonon- 
phonon interaction, whose value is estimated to be 1.8 
cm -1 .— Ho is the phonon self-energy of the doubly de- 
generate Eig mode calculated for n — and ntop = 0, and 
it is given by il = Re [u^n^O, 0)u*] , where k = 1,2 
corresponds to the two E 2g modes and are the corre- 
sponding eigenvectors. In presence of doping charge the 
lattice parameter changes and the corresponding vari- 
ation of the G modes frequencies can be obtained for 
zero electric field by ab initio calculations^ Aw(n) — 
—5.75 10 3 Aa(n) cm -1 , where Aa(n) is the relative vari- 
ation of the lattice parameter, as in Eq.(2) of Ref.— In 
this work we do not include in the phonon calculations 
the effects due to direct interaction of the system with ad- 
sorbates and with the substrate donating doping charge. 
However, as shown for intercalated graphite ) 36 ' 37 we as- 
sume that this is not relevant also for the in-plane vibra- 
tional modes of bilayer graphene. 

The eigenvalues of Q,^ are of the form (u>i + iji/2), 
where u>i is the frequency of phonon i and 7, is the full 
width half maximum (FWHM), given by the EPC and 
the anharmonic phonon-phonon interaction. In the gen- 
eral case of finite n and ntop the four eigenmodes of f2 M „ , 
e % , are still two by two degenerate, but they are a super- 
position of the E2 g and E u eigenmodes of the unbiased 
bilayer graphene. Their relative Raman intensities P R are 
calculated as: 
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where E i / fl= 1 - 

We have shown that the screening properties of bi- 
layer graphene under the application of an external 
electric field are characterized by inter- and intra-layer 
polarizations^ Most of the calculations of the band gap 
as a function of n and n top are based on TB models^^ 
which usually consider only the inter-layer palarization, 
resulting in an overestimation of the band gap. In the 
present work the band structure of the tt electrons [eki 
and |ki) in Eq.fTTJ)] is obtained using the scheme presented 
in Rcfi^ The band gap is computed by ab initio DFT cal- 
culations, including the GW corrections, and both inter- 
and intra-layer polarizations are fully taken into account. 
The full band structure of gated bilayer graphene is then 
computed using a TB model, which is able to reproduce 
all the important features of the DFT calculated, GW 
corrected bands, including the electron-hole asymmetry. 

In order to compute AH^ one needs to calculate the 
derivative of the tight binding Hamiltonian with respect 
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FIG. 2: (Color online) Raman shift in bilayer graphene for 
ntop = ntot and for nt op = ~nt, ot . Calculated values of the 
shifts are connected by lines. For a given value of n (panel a) 
or n t0 p (panel b) there are two phonon modes represented with 
two rectangles. The height of the rectangles is the FWHM 
and the areas are proportional to the relative Raman intensi- 
ties (i.e. the integrated area of each peak) of the two modes. 
Thus, the ratio of the widths of the two rectangles is equal to 
the ratio of the maximum heights of the two Raman peaks. 
When the ratio is less than 0.1, the mode with the smallest in- 
tensity is colored in gray (red) , otherwise is black. Circles are 
experimental results from Ref.— and the errorbars represent 
the experimental FWHM. 



to the atomic positions. However, only the variation of 
the first nearest-neighbors in-plane hopping parameters 
71 turns out to be relevant, and this only term is thus in- 
cluded in AH 11 . The value we use for this quantity is 5.8 
eV A -1 , which derives from the ab initio GW-calculated 
EPC at r for the E^g mode in monolayer graphene^S 



III. RESULTS 

In Fig[2]-a we show the calculated Raman shift as a 
function of n, for the case ntop=nbot, where the external 
electric field and the band-gap are kept fixed to zero, 
as theoretically studied in Ref£L In this case the E^ g 
and E u modes do not mix by symmetry. Analogously 
to monolayer graphene, the Raman active modes show a 
singularity when the Fermi energy is half of the phonon 
energy. In Fig[2]-b we show the calculated Raman shifts 
for the case ntop = —nbot, as a function of ntop- This is a 
special situation where the external electric field is varied 
while n is kept fixed to zero, as realized in recent infra-red 
experiments.—^ In this case, the mixing between E-2, g and 
E u modes is weak due to the antisymmetric allocation of 
charges in the two layers. The Raman active modes show 
a singularity in the frequency and a divergence in the 
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FIG. 3: (Color online) Raman shift in bilayer graphene as a 
function of the electron concentration n, for different values 
of ntop- See the caption of Fig(2] 
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FIG. 4: Ratio between the relative Raman intensities of the 
highest and lowest mode, a, and Raman shift, b, as a func- 
tion of n. Filled dots are experimental results from Ref . 15 ' 33 , 
shifted by n°=- 1.8 10 13 cm . Lines are the theoretical values 
for n top =-0.15, -0.25, -0.35 10 13 cm" 2 . 

FWHM when the band-gap is of the order of the phonon 
energy. 

In the most general situation the electric field and n 
are both finite. In Fig[3] we show the Raman shift of 
bilayer graphene as functions of n, for different values 
of n t0 p- In these cases, the E2 g and E u modes do mix, 
and at certain values of nt op and n two modes become 
Raman visible. Our results show an asymmetry between 
positive and negative values of n top - For instance, in the 
case Tif p 2.4 10 cm the FWHM of the lowest mode 
at n rj -4 10 13 cm -2 is higher than the same quantity 
for the case n top = -2.4 10 13 cm -2 at n ~ 4 10 13 cm -2 . 
This is due to the electron-hole asymmetry in the band 
structure, which is properly described in our calculations. 
Moreover, the asymmetry between positive and negative 
n is enhanced by the effect of the lattice spacing variation 
induced by the doping charge. Our results are qualita- 
tively in agreement with recent calculations^ based on 
TB model. They are however quantitatively different, 
because in our calculations we include the electron-hole 
asymmetry, the lattice spacing variation due to doping 
charge, and both inter- and intra-layer polarizations. The 
dependence on n of the frequencies, FWHM, and relative 
Raman intensities is strongly influenced by nt op , and on 
the basis of this observation we claim that the amount 
of uncontrolled n top and n^t can be estimated from the 
Raman spectra of bilayer graphene when two modes are 
observed. 

We now compare our calculations to the experimental 
Raman spectra of bilayer graphene where the splitting of 
the G mode is reported^ In this work charges are inten- 
tionally induced on the system by applying a gate voltage 
between the bilayer and the Si02 substrate. However, 



unintentional n top and 7i£, ot arising from the atmosphere 
and the substrate can be present at zero gate voltage. 
By comparing the experimental and calculated Raman 
shifts as a function of n for different n top , we estimate a 
total electron concentration at zero gate voltage n° = -1.8 
10 13 cm -2 . In FigdJ-a and -b we show the ratio between 
the relative Raman intensities of the highest and lowest 
mode, and the Raman shifts, respectively, as a function 
of the electron concentration n, for different values of 
n top - The former one strongly depends on n top , while 
the frequency shifts have a weaker dependence. The best 
agreement between theory and experiments indicates an 
unintentional charge coming from the atmosphere n® op =- 
0.25 10 13 cm~ 2 . From our estimate of n°, we deduce an 
unintentional charge from the substrate n^ ot =n°-n^ op = 
-1.55 10 cm . The agreement between experimental 
data and theoretical results is good. However, we notice 
that the slope of the theoretical curves is underestimated 
with respect to the experimental ones. This could be 
possibly due to local desorption of molecules and dop- 
ing variation induced by the laser light, or to hysteresis 
effects in the doping dependence on gate voltage. 

Finally, in the right side of FigJ^-b, we compare 
our theoretical results to the experimental frequencies 
and broadenings from recent infra-red measurements in 
RefJi, where the doping charge is kept fixed to zero and 
the electric field is varied. The agreement is excellent 
with our lower frequency mode. In our calculations, the 
lower mode has a weak projection on E u . However, this 
mode is strongly coupled with the electrons, as testified 
by the large FWHM. Such coupling enhances the effec- 
tive charges associated with E u and increases the infra- 
red activity*i&2i>2£ Indeed, in Fig.3 of Ref^ the mea- 
sured infra-red intensity is maximum when the band-gap 
equals the phonon energy (about 0.2 eV), i.e. when the 
FWHM and thus the coupling of the mode with the elec- 
trons are maximum, while it decreases when the FWHM 
decreases. 



IV. CONCLUSIONS 

In summary, we have computed the Raman spectra of 
gated bilayer graphene, which is strongly influenced by 
the interaction with the environment. We claim that by 
the analysis of the splitting of the G mode in Raman 
measurements it is possible to estimate the amount of 
unintentional charges coming from the atmosphere and 
from the substrate. Here we compare our calculations 
with the only experimental data available on the G mode 
splitting in bilayer graphene, and we give an estimate of 
the unintentional charges coming from the environment 
in this experiment. In order to facilitate the comparison 
of new experimental results with our theoretical calcula- 
tions, we provide as additional material a set of computed 
Raman shifts, FWHM, and relative Raman intensities as 
a function of n, for different values of n top ^ 
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